High energy collisions of strongly deformed nuclei: 
An old idea with a new twist 
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UU collisions can provide about 30% larger densities compared to central PbPb ones. New aspect 
is generation of rather deformed initial states. We show that those can be effectively used to resolve 
a number of outstanding issues, from corrections to hard processes, elliptic flow (the QGP push 
issue), and the mechanism of J/ip suppression. UU collisions are studied by a simple Monte-Carlo 
model, and it is shown how selecting two control parameters - the number of participant nuclei and 
deformation - one can select particular geometry of the collision. 



I. INTRODUCTION 

An "old idea" mentioned in the title is to select the 
head-on "long- long" collisions, by simply triggering on 
maximal number of participants Np. Although it is 
kind of a folklore of the field, the only written mate- 
rial on it I found is a memo ||i| written by P.Braun- 
Munzinger. His estimates show that, due to larger A 
and deformation, the gain in energy density for UU over 
AuAuQcan reach the factor 1.8. Although our study had 
found smaller numbers, new applications are proposed. 
They are mostly related with a different geometry of the 
collisions, the parallel one, which incorporates about the 
same energy density as available in central PbPb colli- 
sions, with significant deformation. 

The general attitude of this study was to look how UU 
collisions can help to understand the existing open issues 
of the SPS heavy ion program. It certainly is of interest 
to RHIC program, although it is probably premature to 
discuss it now. The emphasis is made here on event se- 
lection, and other interesting options (like using targets 
with crystals with naturally aligned U) are not studied. 

Let me start with outlining simple argument for head- 
on collisions. Representing U as a homogeneous ellipsoid 
with one long (Ri) and two short (i?s) semi-axis, one can 
related their ratio to deformation parameter d used in 
nuclear physics (see e.g. B) 



Ri _ 1 + 4(5/3 , 
Rs ^ ^1-2(5/3' 
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For 6u w .27 this ratio is 1.29, the basic deformation 
ratio to be used below. 

It is convenient to think first in terms of "wounded" 
or "participant" nucleons first, a purely geometric con- 
cept, and only then consider real multiplicity (entropy) 
production. (We will follow such logic throughout the 
paper.) Let us thus start with comparing the density 
of participants per transverse area rip = Np/ (ttR^^), for 
"long- long" collisions of the deformed nuclei vs the spher- 
ical one with the same A and R — {R'^RiY^^. The effect 
only comes from reduction of the area, so 
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For A=238 we wiU use Ri = 8.4, R, = 6.5, R = 7.0/m, 
and so deformation alone increase the rip by 1.16. For 
U and Pb (i?p(,=6.78 fm in such model) one gets the 
participant density gain 1.24. 

Transferring this into initial entropy density, one 
should recall that for (spherical) AA collisions 



dN_ 
dy 



(2/ = 0) ~ A 
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with a w .12. So, assuming as usual that final multi- 
plicity is proportional to the initial entropy density, we 
see that there is a correction to the simple idea that each 
participant nuclei gives the same (energy dependent) con- 
tribution to the spectrum. This non-zero a incorporates 
both (i) additional increase in multiplicity, and (ii) ex- 
tra stopping (shift toward mid-rapidity): we are only 
interested in their combination dN/dy(y=0). Further- 
more, it is natural to think that transverse dimensions 
of the system enter trivially here, and so these extra ef- 
fects due to increased density. With this additional factor 
~ rip" we obtain the UU/PbPb total initial density gairw\ 
1.241+3" = 1.34. 

The main questions addressed in this paper are two- 
fold. One is to make some realistic estimates of the effect, 
not just for a particular configuration but for ensemble 
of events selected by some experimentally accessible cri- 
teria. The second is to outline possible applications of 
high energy collisions of the deformed nuclei. 



II. UU COLLISIONS VERSUS PBPB: THE 
SIMULATION 

Simple Monte-Carlo program was written, which ini- 
tializes nucleons inside nuclei and follow their paths 



'Below we use PbPb instead. 



^At fixed time after the crossing, a la Bjorken. If more 
stopping means earlier equilibration, the available density is 
larger. 



through another one. Since we are not really inter- 
ested in peripheral collisions, we did not included dif- 
fuse boundary of nuclei and used the ellipsoids described 
above. We also ignored probabilistic nature of the inter- 
action, considering transverse distance between nuclcons 
R < (cm/Tr)^/^ to be sufficient reason to make both of 
them participants. So, the only source of fluctuationsn 
are random positions of the nucleons inside the nucleus. 

Spherical nuclei have only one parameter - impact pa- 
rameter b - which in such classical treatment determines 
the mean number of participants. Deformed nuclei have 
in general 5 such parameters: b and 4 spherical angles 
9i,(j}i,i = 1,2 indicating the orientation of their longer 
axes at the collision moment. The main objective of the 
calculation is to see how well one can actually fixed those, 
by using experimentally available information. 

Few words about our definitions. After all participant 
nuclei are identified, we calculated the tensor 



J-ij — ^ '^i ^j -^ 
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in transverse plane, diagonalize it and find its eigenval- 
ues i?^ , i?^ . The density of participants we use below 
is defined as n^ = Npart/{T^R+R-) aud deformation as 
R+/R-. 



with corresponding A for for both cases. Triggering on 
large Np (or forward energy) one effectively eliminates 
spectators (and many comphcated geometries). 
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FIG. 2. Distribution of R+/R- (long-to-short semi-axis 
in transverse plane), with selection Np > 0.9(2yl). Solid 
(dashed) are for UU and PbPb collisions 




FIG. 1. Density distribution with selection Np > 0.9(24). 
Solid (dashed) are for UU and PbPb collisions 

If no selection is made, the density of the participants 
Up shown in|l| is very similar, except for the tail region. A 
gain in np of the order of 16% seem indeed possible, with 
reasonable loss of statistics. However, if one introduces 
"centrality" cuts, the situation is different. We would 
define a fair "centrality cut" by restricting Np > 0.9{2A), 



The first striking feature one finds after such cut is that 
distribution of the deformations of the initial 2d ellipsoid 
R+l R^ is very different: see figure Fig.|[ The ratios as 
large as 1.35 are accessible, while in PbPb collisions all 
the "central" collisions are very spherical. The maximal 
deformation, not surprisingly, is of the order of deforma- 
tion of U|^ Those correspond to collisions with two long 
directions parallel to each other and orthogonal to the 
beam. It is about the same as is obtained for medium 
b for spherical nuclei, but at larger energy density and 
larger system (see below). 

The joint distributions in participant density - defor- 
mation plane for (most central) UU and PbPb collisions 
are compared in Fig.H(a,b). The main message one can 
get from it is that strong correlation between the defor- 
mation and density, existing for spherical nuclei, is to 
some extent relaxed for UU. 

How one can measure the 2d deformation of the ini- 
tial conditions? The measured elliptic deformation of 
spectra of secondary particles, pions or nucleons, V2, is 
proportional to this initial deformation (with EOS de- 
pending coefficient) and should have similar distribution, 
divide measured distribution over V2 into more and less 
deformed. 

Suppose now that one uses both control parameters. 



"'"And it is by no means assumed to be accurate account for 
fluctuations. 



^Note that it is comparable to what is obtained for medium b 
for spherical nuclei, but now reached for larger energy density 
and larger system. 



Npart and V2- Is it really possible to select the partic- 
ular geometries of the collisions we want? In Figjj(a) 
one can see that it is to a significant effect correct: the 
less — deformed sample is rich in the region cosO ~ 1, 
or in "head-on" collisions, while the more — deformed 
collisions have none of them, and concentrate at small 
COS0. In Fig.^b) one can see that the same selection 
of events correspond to the difference of the azimuthal 
angles 0i — 02 to be peaked around 0, or have a flat 
distribution, respectively. 
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FIG. 3. Distribution over participant density rip [fm~^] 
vs deformation R+/R^, for (a) UU {Npart > 428) and (b)) 
PbPb {Npart > 374) collisions, respectively. 



FIG. 4. Distributions of original angles for UU collisions, 
Np > 0.9{2A). (a) Distribution in cos{9), the angle between 
the long axis of U and the beam. The solid (dotted) his- 
tograms are for less deformed R+/R- < 1.1 ( more deformed 
R+/R- > 1.2) initial states . (b) Distribution in difference 
between polar angles \(j>i — (/fill, solid (barred) histograms are 
for less and more deformed collisions, same selection. 



Finally note also that the figures presented above have 
shown only the density of participants Up: let us now 



recall the correction factor 



,3a 



we discussed in the In- 



troduction for the initial energy density. It can be seen 
as additional non-linear deformation of the axis, when 
going from Up to dN/dy. The contrast between UU and 
PbPb in dN/dy is larger by another 8%. Note also, that 
we have only discussed local quantities, like participant 
density rip. For many applications the absolute size of 
the system is of similar and sometimes even larger im- 
portance. 



III. POSSIBLE APPLICATIONS 

A. Hard versus soft processes 

The so called hard processes include Drell-Yan dilep- 
ton production, two-jet events, heavy flavor production 
etc. Those are described in the first order by the parton 
model, and so are the simplest to treat for any geome- 
try. For example, for head-on collisions discussed in the 
introduction, the output should be simply 



^hard 



Ityltf 



(5) 



Significant interest is related with various QCD 
corrections to the parton model: some of them are re- 
lated to initial state re-scattering, and some to the final 
state ones. Examples of the former are "shadowing" of 
nuclear structure functions and increase in parton pi, the 
latter processes lead to the so called "jet quenching" or 
energy losses. The debates about their relative contri- 
butions continues: one, at one talk as a semi-joke I pro- 
posed to consider rectangular nuclei with different R^ , Rt 
to separate them . Now we propose a particular realiza- 
tion of this idea. 

How large level arm do we have, with UU collisions? 
For head-on coUisions (discussed above) Rz = Ri.Rt = 
Rs , while for "parallel collisions" with maximal deforma- 
tion it is the other way around, Rz ~ Rs and one of 
the Rt = Ri- So the ratio Rz/Rt varies between about 
1.3 and 1/1.3, or a span of 1.7. This may be enough to 
disentangle different mechanisms. 



B. Ellipticity and EOS 

In high energy collisions the shape of the "initial al- 
mond" for non-central collisions leads to enhanced "in- 
plane" flow, in direction of the impact parameter H . It 
is very important because (as pointed out in |g[(a)) it is 
developed earlier than the radial one, and thus it may 
shed light on Equation of State (EOS) at early time. The 
particular issue is whether we do or do not have QGP at 
such time, at SPS or RHIC. Recent review of ellipticity 
one can find in [Q. 

Ellipticity is now measured by the asymmetry of the 
particle number, or Vi harmonics defined as 



—- = — + — cos{2(j)) + — cos(4(/i) 

ad) ZTT TT TT 
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rather than asymmetry of the momentum distribution, 
which was used more at low energies. Furthermore, Vi 
are often additionally normalized to the spatial asymme- 
try of the initial state (the "almond") at the same b, 
"2 = (-R+ — R-)/iR% + R-) in order to cancel out this 
kinematic factor and to see the response to asymmetry. 
There are two ways to look at elliptic flow: by using 
collision energy dependence, or centrality dependence at 



fixed beam. The first method is more difficult, but it 
deals with fixed geometry and only slowly changing size 
of the system and relevant densities. For recent discus- 
sion of it see ref.3. The centrality dependence is easier 
to measure, but interpretation is more complicated, be- 
cause increasing b we make "almond" more elliptic, but 
also much smaller and thinner: eventually finite size cor- 
rections reduce pressure build-up. This is of crucial im- 
portance for the "QGP push" issue at SPS, since it should 
be present only for the largest densities available, and at 
small b the V2 is very small and difficult to measure. 

Preliminary NA49 data presented at QM99 g] may 
indeed indicate "the plasma push" , as enhancement of V2 
at small b. However much more work is clearly needed 
to understand this complex interplay of EOS and finite 
size effects. 

U collisions discussed below in principle provide the 
means to decouple finite size and deformation issues. In 
particularly, the deformation of U (about 1.3) is actually 
enough to generate well measurable V2 of the order of 
several percent, without significant loss in density. 



C. J/tp suppression 

One of the most intriguing observables, the J ftp sup- 
pression, is unique in its significant centrality and A- 
dependence. Here is not a place to discuss it in any 
details, but let me make few remarks about recent de- 
velopments. 

New NA50 data reported at QM99 have clarified the 
situation for the most central collisions: using now very 
thin target, it was found that 1996 data suffered from 
multiple interactions. In fact there is a significantly 
stronger suppression at small b. Furthermore, it seem 
like the two component picture, with separate x ^nd ip 
thresholds, really emerges. In view of this, it is desirable 
to increasing the density and/or the famous variable L: 
only deformed U provides an opportunity here at SPS. 

In order to discriminate experimentally different ideas 
on the nature of J/ip suppression, we should be able 
to tell whether suppression happened quickly or need a 
longer time. The old idea is to study suppression depen- 
dence on pt- Unfortunately changing pt we also change 
the kinematics: e.g. destruction by gluons or hadrons go 
better if pt grows. 

Maybe better idea Q is to use azimuthal dependence of 
the suppression. Instantaneous suppression should show 
no asymmetry, but if it takes few fm/c the anisotropy 
should show up. However, as for the elliptic flow, the 
problem is the initial "almond" at b< 8 fm is not very 
anisotropic, and for larger b there is no anomalous sup- 
pression. 

Here too the deformed U can help, providing (with trig- 
ger conditions discussed above) a variety of geometries, 
including the "parallel collisions" . 
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